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Abstract 



^^ ' One of the fundamental principles driving diversity or homogeneity in domains such as cul- 



tural differentiation, political affiliation, and product adoption is the tension between two forces: 
influence (the tendency of people to become similar to others they interact with) and selection 
(the tendency to be affected most by the behavior of others who are already similar) . Influence 
tends to promote homogeneity within a society, while selection frequently causes fragmentation. 
OO ■ When both forces are in effect simultaneously, it becomes an interesting question to analyze 

^^ I which societal outcomes should be expected. 

In order to study the joint effects of these forces more formally, we analyze a natural model 

^H I built upon active lines of work in political opinion formation, cultural diversity, and language 

r ^ ■ evolution. Our model posits an arbitrary graph structure describing which "types" of people can 

• I influence one another: this captures effects based on the fact that people are only influenced by 

(J • sufficiently similar interaction partners. In a generalization of the model, we introduce another 

graph structure describing which types of people even so much as come in contact with each 

other. These restrictions on interaction patterns can significantly alter the dynamics of the 

^ I process at the population level. 

00 ' For the basic version of the model, in which all individuals come in contact with all others, 

^O I we achieve an essentially complete characterization of (stable) equilibrium outcomes and prove 

\l ■ convergence from all starting states. For the other extreme case, in which individuals only 

come in contact with others who have the potential to influence them, the underlying process is 
g : signiflcantly more comphcated; nevertheless we present an analysis for certain graph structures. 
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1 Introduction 

Selection and Influence. Human societies exhibit many forms of cultural diversity — in the 
languages that are spoken, in the opinions and values that are held, and in many other dimensions. 
An active body of research in the mathematical social sciences has developed models for reasoning 
about the origins of this diversity, and about how it evolves over time. 

One of the fundamental principles driving cultural diversity is the tension between two forces: 
influence and selection. Influence refers to the tendency of people to become similar to those with 
whom they interact, whereas Selection is the tendency of people to interact with those who are 
more similar to them, and/or to be more receptive to influence from those who are similar. 

Both of these forces lead toward outcomes in which people end up interacting with others like 
themselves, but in different ways: influence tends to promote homogeneity, as people shift their 
behaviors to become alike, while selection tends to promote fragmentation, in which a society 
can split into multiple groups that have less and less interaction with each other. Research that 
offers qualitative analyses for issues such as consensus-building, political polarization, or social 
stratification can often be interpreted through the lens of this influence-selection trade-off. [9l [121 114)^1 
The trade-off between influence and selection, and the development of data analysis techniques 
to try separating the effects of the two, have been integral to understanding and promoting the 
adoption of products and behaviors in social networks [21 [71 [131 [18] , ^^ active line of work at the 
interface of computing, economics, and statistics. 

When both influence and selection are operating at the same time, how should we reason about 
their combined effects? Several lines of modeling work have approached this question, all starting 
from similar underlying motivations, but developing different mathematical formalisms. 

1. Research on political opinions has studied populations in which each person holds an opinion. 
The opinion is represented by a number drawn from a bounded interval on the real line M^, 
or from a discrete set of points in an interval. (For example, the interval may represent the 
political spectrum from liberal to conservative.) Each person is influenced by the opinions 
of others who are sufficiently nearby on the interval, thus capturing the interplay between 
influence (people are shifting their opinions based on the opinions of others) and selection 
(people only pay attention to others whose opinions are sufficiently close) [6l [lOl [11] . 

2. Axelrod proposed a model of cultural diversity in which there are several dimensions of 
culture, and each person has a value associated with each dimension (e.g., a choice of language, 
religion, or political affiliation). Agents are more likely to interact when they agree on more 
dimensions; when two people interact, one person randomly chooses a dimension in which 
they differ, and changes his value so that they now match in this dimension [3]. For example, 
two people who have passions for similar sports and styles of food may end up having an easier 
time (and more opportunity for) associating, and hence an easier time influencing one another 
along another dimension such as religious beliefs. Again, the model represents an influence 
process in which the interactions are governed by selection based on (cultural) similarity. 
Axelrod's model has generated a large amount of subsequent work; see [8] for a survey. 

3. Finally, Abrams and Strogatz exhibited some of the interesting effects that can occur even 
when there are only two types of people. They modeled a scenario in which people speak 
one of two languages. People mainly interact with speakers of their own language, but there 



^The term homophily is often used to refer to the mechanism of selection 14 . However, in other contexts, it is 
used to refer to the broader fact that people tend to be similar to their neighbors in a social network, regardless of 
the mechanism leading to this similarity. Hence, we use the more specific term selection here. 



is gradual "leakage" over time as speakers of one language may convert to become speakers 
of the other [1]. The Abrams-Strogatz model has also generated an active line of follow-up 
results, including explorations of its microfoundations through agent-based simulation [20j 
and analyses of the spatial effects and population density [16]. 



Commonalities among Models. Although the models described above differ in many details, 
they have the same underlying structure: the population is divided into a set of types (the opinions, 
the cultural choices, the language spoken), and a person of any given type may be influenced to 
switch types, but only by others whose types are sufficiently similar. (In the case of the Abrams- 
Strogatz model, there is a preference for one's own type, but since there are only two types, all 
types can influence each other.) This process generates a "flow" as people migrate among different 
types, and we can ask questions about both dynamics (which outcomes the process will reach) and 
equilibria (which outcomes are self-sustaining, in the sense that the flows between types preserve 
the fraction of people who belong to each type). Following the language around Axelrod's work, 
we will refer to this type of process as representing the cultural dynamics of the population. 

In addition to their similarities in structure, these cultural dynamics models also agree in their 
broad conclusions. In the first two models, the population gradually separates into distinct "islands" 
in the space of possible types; subsequently, no further interaction between the islands is possible. 
In the Abrams-Strogatz model, with just two types, the only outcomes that are stable under 
perturbations are the two extreme outcomes in which everyone ends up belonging to the same 
type. Typically, there is also an unstable equilibrium in which each language is spoken by a non- 
zero fraction of the population. 

The most salient difference among the models is the structure that is imposed on the set of 
types. In each case, there is an undirected influence graph T on the set of types: when a person 
of type u interacts with a person of type v, the person of type u has the potential to switch to 
(or move towards) v provided that u and v are neighbors in T (i.e., provided that u and v are 
sufficiently similar according to the interpretation of the model). In the models of one-dimensional 
opinion dynamics on a discrete set, the graph T is the k power of a path for some k > 1 (types 
are similar enough when they are within k steps on the path); in Axelrod's model, the graph T 
is the k power of a (not necessarily binary) hypercube. The Abrams-Strogatz model shows that 
these kinds of processes can exhibit subtle behavior even on a two-node influence graph T. 

The present work: Cultural Dynamics on an Arbitrary Influence Graph 

All of the prior results apply only to highly structured, symmetric graphs (essentially hypercubes 
and paths), whereas in some of the settings that the models seek to capture, the set of types can 
have a less orderly structure. A basic open problem is to analyze cultural dynamics on an arbitrary 
influence graph. 

This is the problem we address in the present work, where we develop techniques for resolving 
some of the main questions on arbitrary graphs. For a natural formulation of cultural dynamics 
on an arbitrary influence graph (which refer to as the global model, for reasons explained later), 
we prove convergence results and precisely characterize the set of all stable equilibria. We then 
consider generalizations of the global model, extending some of our convergence and stability results 
to these more general settings and posing several open questions. 

The Global Model. We now describe the global model in more detail. Because the models from 
the earlier lines of work discussed above differ in many of their details, there is no meaningful way to 
simultaneously generalize all of them in a precise syntactic sense. Instead, our goal is to formulate 



a version of cultural dynamics that exhibits the same basic interplay of selection and influence — 
specifically, the idea that influence only happens among types that are "close together" — while 
allowing for an arbitrary graph on the set of types. 

Let T be a graph on a set of types V; for each type u £ V, let T„ C y denote the set of 
ti's neighbors in T. As is standard in many of the approaches to cultural dynamics, we model 
the population as a continuum: at the start of the process, each type u a V has a non-negative 
population mass associated with it, corresponding to the fraction of the population that initially 
has this type. (Consider, for example, the fraction of the world's population that belongs to a 
certain religion or speaks a certain language.) Time evolves in discrete steps t = 0, 1,2, . . ., and 
Xu{t) denotes the mass on type u at time t. The full state of the population at time t is thus given 
by the mass vector x{t), the vector of values Xu(t) for all u G V. 

We define a discrete-time dynamical system in which x„(t + 1) is determined in terms of the 
mass vector x{t). The dynamical system is motivated by imagining that each person chooses a 
random other person to interact with. Selection effects are captured in two ways by the model: 
first, people are more likely to interact with their type; and second, they only have the potential to 
be influenced when they interact with an individual of their own or a neighboring type. Specifically, 
each person is a times more likely to choose an interaction partner of their own type than someone 
of a different type, for a parameter a > 1. When a person of type u chooses to interact with a 
person of type v, such that v G T^, then with probability p, he will switch to type v, where p € (0, 1] 
is a fixed parameter. To express this dynamic numerically, we let Muit) = ci^uit) + Yly^Y\Su} Xvit)- 
A person of type u chooses to interact with his own type with probability aXu{t)/A4u{t), and with 
any other type v j^ u with probability Xv{t)/Aiuit)- Since a person can change his type from u 
to V only when v £ Tu, the probability that a person of type u keeps his current type after one 

mter action is [1 — p) + p ■ . , , , . 

Mu{t) 

It is easy to translate this idea into a deterministic update rule on the continuum of people: 
following the motivation above, we set 

.„(i + 1) ^ (1 - rt .„(*) + P-^t) . °'^°''' + ^mm"""' "'"' + p E ^'^"(*) ■ lif)' w 

where the first two terms represent the mass that stays at u from time t to t+1 (because this fraction 
of u was not influenced or interacted either with other members of u, or with types v ^ {u} U T„), 
and the third term incorporates all the mass that moves from other types v to u. This defines the 
update rule for the mass vectors x in the dynamical system. 

It is natural to think of the parameter p as generally being very small, since most interactions 
between people do not lead to a change of type. However, as it turns out, the value of p does not 
have a major qualitative effect on our results. This is not surprising, since introducing p < 1 (as 
opposed to p = 1) just slows down the "flow" between any two types by a factor of 1/p. We include 
p in the model in order to capture the range of possible speeds at which transitions can happen. 
For example, if the types in our model correspond to dialects of a language, we can choose a small 
p (since the probability that a person changes his dialect is very small). However, if instead the 
types represent opinions in the period before an election, people may switch much more rapidly, 
and a larger p is appropriate. 

Convergence, Equilibria, and Stability in the Global Model. Our first result is that for 
any influence graph T and any initial mass vector x the system converges to a limit mass vector x* . 
We prove this by establishing a system of invariants on the population masses over time, capturing 



a certain "rich-get-richer" property of the process — essentially, that the types of large mass will 
tend to grow at the expense of the types of small mass. 

We next consider the equilibria of this model: we say that a mass vector x is an equilibrium 
if it remains unchanged after one application of the update rule. It is easy to construct examples 
of equilibria that are not stable, in the sense that an arbitrarily small perturbation of the masses 
X* can — after further applications of the update rule — push the masses far away from the 
equilibrium. Such equilibria are less natural as predicted outcomes of the cultural dynamics being 
modeled, since the population would be unlikely to hold its position near this equilibrium. 

To make this statement precise, we use the notion of Lyapunov stability. We say that an 
equilibrium x* is Lyapunov stable (or, more briefly, stable) if given any e > 0, there exists a. 6 > 
such that if ||x(to) — x*\\ < 6, then \\x{t) — x*\\ < e, for all t > ioO For simplicity, we use the Li 
norm instead of the L2 norm in what follows, although this is not an important distinction. 

We prove that x* is a stable equilibrium if and only if the set of active types A{x* ) = {n : x* > 0} 
is an independent set in the influence graph T. The proof is based on the rich-get-richer properties 
of the process; these properties are used to show that after a sufficiently small perturbation to the 
population masses, the amount by which any type with positive mass can grow is bounded. 

Interpretations of the Basic Results. These results establish, first of all, a precise sense in 
which the natural equilibria the model tend to break the population into non-interacting islands. 

In addition to offering a qualitative statement about fragmentation of opinions, the results also 
suggest a way of reasoning about the phenomenon by which opinions on different issues tend to 
become aligned, with an individual's views on one issue providing evidence for his or her views on 
another J17tll9j. To take a concrete example that already appears on the 2-dimensional hypercube 
(i.e. the 4- node cycle), consider a setting in which each individual has either a liberal or conservative 
view on fiscal issues and either a liberal or conservative view on social issues. If we assume that 
people only influence each other when they agree on at least one of these two categories of issues, 
then the graph on the set of types is a 4-node cycle. Since our results on stable equilibria indicate 
that independent sets are favored as outcomes, we can interpret the conclusion in this example as 
predicting that either the whole population will converge on a single node (representing a uniform 
choice of views), or on an independent set of two nodes, in which case an individual's opinion on 
fiscal issues has become correlated with his or her opinion on social issues. 

It is also instructive to compare our results to the main result of pj discussed above. Recall that 
they consider the influence graph T = K2 (two connected nodes) , and they flnd that the two stable 
equilibria are the outcomes in which all the population mass is gathered at a single node. The 
family of dynamical systems they consider strictly subsumes ours in the special case of a two-node 
graph, but for the speciflc system we study, our results imply that their basic finding extends to 
arbitrary graphs: in any graph, the stable equilibria correspond to the non-empty independent sets, 
just as Abrams and Strogatz showed for the two-node graph K2. 

A Generalization: Limiting both Interaction and Influence. We now discuss a natural 
generalization of the model that is significantly more challenging to analyze. In the global model, 
the members of type u can interact with members of all other types, even though they are infiuenced 
only by the types in T^. However, there are settings in which it is more natural to assume that the 



^ One could ask about stronger notions of stability, in particular, asymptotic stability, which requires that there 
exists a 5i > such that if l|a:(to) — a;*|| < 5i, then x{t) -^ x* as t ^ cxd. Asymptotic stability is not a useful 
definition for our purposes; for example, if the underlying influence graph T has no edges, then any assignment of 
population masses is an equilibrium, but none are asymptotically stable, since there is no way for a small perturbation 
to converge back to the original state. On the other hand, all equilibria are Lyapunov stable in this simple example. 



members of a type only ever interact with members of a subset of the other types. For example, 
this may be a reasonable assumption when types represent different languages. 

To capture this idea, we now assume that there are two potentially distinct graphs on the set of 
types V: the influence graph T (as before), as well as an undirected interaction graph S, where T is 
a subgraph of S. Rather than interacting with a person chosen from the full population, a member 
of type u selects an interaction partner from the set 5„ of tt's neighbors in S. It is straightforward 
to write the new update rule for this more general dynamical system, by summing over types in 
Su instead of F \ {u}. Specifically, we can now define A4u{t) = axu{t) + J2vgs ^v{t), and get the 
following generalized update rule: 

Xu[t + I) = {I - p)xu{t) + p-Xu{t) -—J- +p 2^ x^{t)- (2) 

The global model is simply the special case in which the interaction graph S is the complete graph. 
The name global model emphasizes that each type interacts "globally," with all other typeso 

The behavior of this general model is significantly more complex than the behavior of the 
global model; for instance, for arbitrary S and T, it is not even clear whether the process will 
always converge. Intuitively, much of the difficulty comes from the fact that when we consider 
two neighboring types u and v, the sets of types that they are interacting with, 5,^ and Sv, can 
be quite different, whereas in the global model they are both the full set V. Among other things, 
this can lead to violations of the rich-get-richer property that was so useful for reasoning about the 
dynamics of the global model. 

For the general model, we first establish a necessary condition for equilibria, as well as sufficient 
conditions for convergence and stability. We then focus further on the special case in which S = T- 
This is in a sense the opposite extreme from the global model; instead of making S as large as it 
can be, we make it as small as possible subject to the constraint that it contains T as a subgraph. 
Accordingly, we refer to the case S = T as the local model. There are many interesting open 
questions surrounding the behavior of the local model; we make progress on these through initial 
convergence results and the identification of a large class of equilibria that are stable for all a > 1: 
non-empty independent sets for which all nodes in the set are at a mutual distance of at least three. 
In fact, this is an "if and only if" characterization for an important class of influence graphs: those 
whose connected components are trees or, more generally, bipartite graphs. 

An interesting observation is that the local and global models can have genuinely different 
behaviors starting from the same initial conditions: Figure [T] shows an example of an initial mass 
distribution on the three-node path for which the global model converges to an outcome in which 
the mass is divided evenly between the two endpoints, while the local model converges to the 
outcome in which all the mass is on the middle node. 

At a higher level, formalizing the distinction between interaction (5) and influence (T) is a 
potentially promising activity more broadly, particularly in light of the considerable recent interest 
in the effects of information filtering on the political process. (See |15^ [2T] for popular media 
accounts, and [5J for recent experimental research.) The concern expressed in all these lines of 
work is that personalization on the Internet makes it possible to sharply restrict the diversity of 
information one sees, and thus risks accentuating the degree of polarization and fragmentation in 



^ There are clearly many other potential generalizations which could incorporate notions of non-uniform inter- 
action, including different interaction strengths between different pairs of types. Such extensions would lead to 
interesting questions as well, fn the present work, we focus on the generalization with two unweighted graphs 5 
and T because it captures in a direct way some of the additional complexity that is introduced by simultaneously 
modeling limited interaction and influence. 



t>HihHi) (iK^®^i) ®^S^ 

(a) ImUal fractions. (b) Global model (a = 2). (c) ^oca/ mode/ (a = 2;. 

Figure 1: An instance in which different models predict convergence to different equilibria. The global 
model predicts an outcome in which two non-interacting types survive (polarization), whereas the local 
model predicts that only a single types survives (consensus). 

political discourse — essentially, the risk is that people will only ever be exposed to those who 
already agree with them, making any kind of consensus almost impossible to achieve. 

In this context, our general model also brings into the discussion the interesting contrast between 
interaction and influence. Personal filtering of information by Internet users can restrict the set 
of people they interact with (affecting the sets Su), and it can also, separately, restrict the set of 
people who may be able to influence them (affecting the sets T„). These two different effects are 
often bundled together in discussions of information filtering; it will be interesting to see whether 
treating them as genuinely distinct can shed additional light on this set of issues. 

2 Observations on the General Model 

In this section, we develop several observations that apply to the fully general model with an 
arbitrary interaction graph S. In the subsequent sections, we utilize these observations to analyze 
the global model (where S is the complete graph) and the local model (where <S = T). 

We begin by restating the network dynamics in terms of flows on edges. Recall from the 
introduction that the fraction of people of type u who switch to type u € T^ is p-Xv(t)/A4u{t), 
where 7W«(t) = axu{t) + YlveS ^v{t)- Thus, the fraction of the entire population which is moving 
from n to t; is p-Xu{t)x.u{t) / M.u{t)- At the same time as this mass oi p-Xu{t)xy{t) / Ai^^^it) is moving 
from uto v^ a, mass oi p-Xy{t)xu{t) / Mvi^) is moving from v to u. These movements partially cancel 
each other out, and motivate the following definition of the (directed) flow on the edge {v^u) G T: 

h^uit) = pMt) Xuit) (;^ - ;^) . (3) 

The updated mass can then be written as 

Xu{t + 1) = Xu{t) + Y^ fv^u{t). (4) 

It is easy to check that ([4]) is equivalent to the original dynamics ([2]). Interestingly, the M.u{t) 
terms emerge as crucial quantities that determine the direction of the flow; for that reason, we will 
call A4u{t) the interaction mass of node u at time t. 

We say that a node u is active at time t if Xu{t) > 0, and inactive if Xu{t) = 0. We occasionally 
refer to a node u as x-active if it is active in x and x-inactive otherwise. The set of all active nodes 
under x is denoted by A{x). Much of our analysis concerns the structure of the subgraph Tacti^) 
of the influence graph T induced by the active nodes A(x). We begin by characterizing when mass 
vectors are in equilibrium. 

Proposition 2.1 A vector x* is an equilibrium if and only if each connected component C of 
Tactix*) has the property that all nodes u ^ C have the same interaction mass A4u- 



Proof: For the if direction, let (u, v) € T be an arbitrary edge. If (u, v) lies inside a component C 
of 71ct(a^*)) then by assumption, both u and v have the same interaction mass, so the flow between 
them will be 0. Otherwise, at least one of it, v is inactive, so that XuXv = 0, and again the flow is 
0. Thus, X* is an equilibrium. 

For the converse direction, assume that for some connected component C of 7Ict(a^*)) the 
nodes do not all have the same interaction mass. Let C/ C C be the set of all nodes u G C/ 
minimizing Mu{t)- By assumption, U C- C, and because C is connected, there must be a node pair 
u & U,v & C \U which is adjacent in T. Fix this node pair. 

By definition of U, the flow along any edge (u, v) G T is non-negative, and the flow along (u, v) 
is strictly positive. Thus, after one application of the update rule, the mass on node u strictly 
decreases, which means that x* cannot be an equilibrium. D 

The following useful lemma relates convergence and the change in directions of flows: 

Lemma 2.2 If there exists a time to such that the flows do not change direction after time t^, then 
the system converges. 

Proof: Let G be the directed graph obtained by directing each edge (u, v) of T according to the 
direction of the corresponding flow fu-^^ito). By the assumption, these directions stay constant 
after time to- As flow always goes from types with smaller interaction mass to types with larger 
interaction mass, G must be acyclic. Let wi, t;2, . . . , fn be a topological sorting of the graph, so that 
all directed edges of G are of the form (vi,Vj),i < j. 

We define Xi^{t) = J2i=i ^Viit) to be the total mass at time t > to on the k first nodes in the 
topological sorting. Because the total mass in the system is constant, and all flow goes from nodes 
with lower indices to nodes with higher indices, each of the Xi.(t) must be non- increasing as a 
function of t. Since they are also lower-bounded by 0, each Xk{t) must converge to some value Zk 
as t — )■ oo. Therefore, each x„^(t) converges to Zi — Zj_i as t — > cx). D 

Recall that we are interested in characterizing stable equilibria. We next provide a sufficient 
condition. 

Proposition 2.3 An equilibrium x* is stable if it satisfies the following two properties: 

1. The active nodes form an independent set in the influence graph T ■ 

2. The interaction mass of every active node is strictly greater than the interaction mass of each 
of its inactive neighbors in the influence graph T. 

Proof: Let x* be an equilibrium for which both properties hold. A node u is called x* -active if 
it is active in x* and x* -inactive otherwise. Let A be the set of all x*-active nodes, and A^* denote 
the interaction mass of node u with respect to x* . Define 

5 = ^ min (Ml- Ml) > (5) 

uGA, vfA, {u,v)£T 

by the second property of x*. To show stability, we prove that whenever ||x* — x(to)||i < S, the 
system will satisfy ||x* — x(t)||i < 6 for all time steps t > to- 

The key step of the proof is to establish that for each node u (z A, the mass x„(t) is non- 
decreasing over time. For contradiction, assume that t is minimal, and u G A an x*-active node, 
such that Xu{t + 1) < x„(t). Because A is an independent set in T, all of li's neighbors in T must 
be x*-inactive. By definition of t, we have that "^.^ffj^Xvit) < 6, and ^^,g^2;t)(t) > ^j,g^ a;* — 6. 



In particular, this implies that A4u{t) > A4* — a5 and for any neighbor v of u, ^Ay{t) < M*^ + a5, 

so we obtain that Mu{t) - M^{t) > {Ml - Ml) - 2a6 > (2a + 1)6 - {2a)6 > 0. 

This implies that no flow could have been directed from u to v, for any w € T„, contradicting 
that w's mass decreased. 

Finally, because each of the Xu{t),u € A is non-decreasing, mass can only move among x*- 
inactive nodes, or from x*-inactive nodes to x*-active ones. Therefore, it holds that J2uga \^u{t) — 
Xu{to)\ = T.uiAXu{to) - T^u^A^uit). Thus, 

\\x{t) -X*\\i < ^ \Xu{t) -Xuito)\ + ^ \Xu{to) - xl\ + y^Xujt) 
ueA u£A u^A 

= '^\xu{to) - xl\ + '^Xu{to) = \\x{to) -X*\\i < 6, 
u€A u^A 

SO the system is stable. D 

3 The Global Model 

In this section, we analyze the global model. The definition of the general model states that flows 
are always directed from nodes with smaller interaction mass to nodes with larger interaction mass. 
For the global model, this property is simplified significantly: flow is always directed from types 
with smaller mass to types with larger mass. This property lets us achieve an almost complete 
understanding of the global model. We show that for this model, the system always converges, 
and we present a complete characterization of which equilibria are stable. First, we characterize 
equilibria by applying Proposition 12.11 to the global model. 

Corollary 3.1 Under the global model with a > 1, the system is at equilibrium, x* if and only if 
the following holds: for every connected component C of T&ct{x*), all nodes u G C have the same 
mass. 

Proof: Proposition 12.11 guarantees that x* is at equilibrium if and only if for each edge (n, n') G 
Ta,ct{x*): axl + X]„g5 xl = axl, + J2vgS , ^v- ^^ the global model, for any node u, the set Su 
consists of of all types but u itself, implying that the sum cancels out, and we obtain (a — l)x* = 
(a — 1)3;*/. For a > 1, this implies x* = x*/. D 

We next show that the system always converges; the proof relies on the key invariant that for 
any 1 < k < n, the total mass of the k smallest types never increases over time. More formally, we 
define the following quantities: 

Definition 3.2 Let yi{t) < y2{t) < . . . < yn{t) be the node masses sorted in non- decreasing order. 
Define 

Yk{t) = y^yi{t) = min Vx„(t) (6) 

to be the sum of the masses of the k smallest nodes at time t. 

The following lemma formally captures the notion that the rich get richer in the global model. 
Lemma 3.3 For every k, the function Yk{t) is non-increasing in t. 



Proof: Let t,k be arbitrary, and assume for contradiction that Yk{t + 1) > Yk{t). Let R{t) be 
the set of k nodes achieving the minimum in (j6]) at time t, and similarly for R{t + 1) at time t + 1 
(breaking ties arbitrarily). By assumption and the definition of R{t + 1), we have 

Y^ x^{t) < Yl xvit+i) < Yl ^^(*+i)' 

v€R{t) veR{t+l) veR{t) 

so the total mass of nodes in R{t) must have strictly increased from time t to time i + 1. In 
particular, because the total mass of all nodes is constant over time, this means that the flow 
from R{t) (the complement of R{t)) to R{t) was strictly positive, so it must have been strictly 
positive along some edge {u,v) with u ^ R{t),v € R{t). But by definition of R{t), we have that 
Xu{t) ^ Xy{t), a contradiction to the flow dynamics of the global model. D 

Theorem 3.4 Under the global model, the system converges for any influence graph and any start- 
ing mass vector x(0). 

Proof: By Lemma [3.31 each function Yj{t) is non-increasing in t. As all masses are non-negative, 
the Yj{t) are also bounded below by 0. Hence, each function Yj{t) must converge to some value 
Zj. Thus, each function yj{t) must converge to Zj — Zj^i =: Zj. It remains to show that this also 
implies convergence of x{t). 

Let (5 > be at most the smallest difference between any two distinct Zj, i.e., 5 < minjj-.^.^^. \zi — 

Zj\. 

Let to be large enough that \yi{t) — Zi\ < ^ for all t > t^. Note that this also implies that 
y3it)-Z,<T^,=rMt)-Zi\<j-4^. 

We will show that if for some node f , there exists a time t > to such that Xv{t) = yj{t) and 
Xv{t + 1) = yj'(t + l), then Zj = z'y Assume towards a contradiction that this is not the case; hence, 
there exist v and t> to such that Xt,(t) = y^it)^ Xyit + 1) = yj'{t + 1) and Zj ^ Zj. Furthermore 
assume that j < j'. (The proof for the case of j' < j is similar.) Then, Xy{t) < Zj -\- -^ and 
x^{t + 1) > Zj' — 2^ > -Zj -| — §^^- Therefore, in time step t, v gained at least a mass of ^~ 6. 

The source of this mass can only be nodes which at time t had smaller masses than v. Hence, 

2r7 — 1 () 2n — ^ 

y,_i(t + l)-Z,_i < y^._i(t)-z,_i-^— 5 < ij-i).--^—6 < 0. 
This is a contradiction to the fact that l^_i(t) is non increasing and converges to Zj^i. D 

3.1 Characterization of Lyapunov Stable Equilibria 

For the global model, the properties required for Proposition 12.31 hold for any independent set, since 
the interaction mass of active types is always greater than the interaction mass of inactive types. 
Therefore, any equilibrium in which the set of active nodes is independent is stable. To complete 
the characterization, we show that the converse is also true. 

Theorem 3.5 In the global model, an equilibrium x* is Lyapunov stable if and only if the active 
nodes form an independent set. 

Proof: It remains to prove the "only if" direction. Assume the active nodes in an equilibrium 
x* do not form an independent set. We will prove that x* is not Lyapunov-stable. 

Let C be a connected component of size \C\ > 2 in Tact{x*)- By the assumption that the active 
nodes in x* do not form an independent set, such a connected component exists. Notice that each 
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component of Tg^ctix*) evolves in isolation, so we can focus on only C for the rest of the proof. 
Therefore, by Corollary 13.11 x* = // for all v € C, for some value fi. 

Let u,v € C be two arbitrary nodes, and 6 > he arbitrarily small. Consider the following 
perturbation: x^j = x* + 5, x„ = x* — 6, and x^ = x^ for all w ^ u,v. By Theorem 13.41 the system, 
starting from the perturbed vector x, will converge to some new equilibrium y. By Lemma 13.31 
the smallest mass of any node in C will always be at most n — 6 during the process. All y-active 
nodes must have the same mass; therefore, if all nodes were active in y, they would all have to have 
mass at most fi — 5, which would imply that mass has disappeared from C, a contradiction. Hence, 
at least one node of C must end up inactive in y. In particular, this means that ||x* — y\\ is not 
bounded in terms of 6, and x* is not Lyapunov stable. D 

4 The Local Model 

In the previous section, we have given essentially complete characterizations of convergence and 
stability of equilibria under the global model, in which all types have the potential to interact, even 
though only certain pairs of types can influence each other (according to the graph T). 

We now consider the local model, which is at the other extreme of our general family: here, the 
interaction graph S is the same as the influence graph T; hence, interactions occur only between 
individuals who also have the potential to influence each other. (We will generally denote this 
underlying graph by T, with the understanding that S = T.) We find that the problems of 
convergence and stability are much more challenging in this case. For the global model, we were 
able to extract very useful organizing structures in the dynamical system that gave us a natural 
progress measure toward convergence. But as is well known, in general, a dynamical system on 
even a small number of variables may have convergence properties that are extremely difficult to 
analyze or express. Given the complex behavior of the update rules for the local model, we find 
that the convergence and stability questions are already difficult on graphs T with a small number 
of nodes, and we focus our results here on such cases. (Of course, based on the motivating premise 
of the model, even systems with a small number of variables are frequently natural, corresponding 
to selection and infiuence dynamics in societies with, for example, a small number of languages, a 
small number of political parties, or a small number of dominant religions or cultures.) 

We begin by considering the case q > 1 and first prove the following theorem: 

Theorem 4.1 Under the local model, if the influence graph is a 3-path, then the system converges 
from any starting state. 

The full proof is provided in Appendix [Aj but we provide a brief outline here. The subtle 
difficulty arises due to the fact that the flow between two types u and v does not necessarily go 
in the same direction over all time steps, but instead may change its direction. To keep track of 
the changes in direction, we define a configuration of the system to be a labeling of all edges (n, v) 
in T by the direction along which flow is traveling (i.e., whether from u to v or from v to u). In 
the case of a 3-node path, there are four possible configurations. We study transitions among the 
configurations as the system evolves over time; we show that each configuration is either a sink, 
which cannot transition to any other configuration, or it has the property that any change in the 
direction of an edge leads to a sink configuration. This ensures that there can be at most one change 
in the direction of fiow as the system evolves; hence, there is a time to such that for any t > to, no 
flow changes its direction. After this point. Lemma 12.21 guarantees that the system converges. For 
the case a > 2, we show this fact only for the 3-path; for a < 2, we establish a more general result, 
showing the same fact for arbitrary star graphs. 
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For a = 1, we are able to prove convergence if the active subgraph is a path of n < 5 nodes. The 
proof requires different techniques than the ones we use for a > 1: for paths of more than 3 nodes, 
flows on edges can change their direction infinitely often. The proof is provided in Appendix [Bl 

4.1 Characterization of Universally Stable Equilibria 

Next, we turn our attention to Lyapunov-stable equilibria. We focus on a very strong notion of 
stability: stability for all a > 1. Formally, we call a mass vector x a universally stable equilibrium if 
a; is a Lyapunov-stable equilibrium for every a > 1. Our goal here is to investigate which equilibria 
are universally stable. Such equilibria are robust to (a very idealized notion of) a changing the 
environment, as expressed by varying a. 

Our main result for universally stable equilibria is a complete characterization for influence 
graphs that are forests, and more generally for influence graphs whose connected components are 
bipartite graphs. 

Theorem 4.2 Consider the local model with a > 1. AssuTne that all connected components of the 
influence graph T are bipartite graphs. Then, a mass vector x* is a universally stable equilibrium 
if and only if the distance between any two active nodes in x* is at least 3. 

The proof of Theorem 14 . 21 consists of several sub-results, all of which hold for arbitrary influence 
graphs, and imply the desired characterization under the assumption in the theorem. It is worth 
noting that these sub-results constitute significant progress towards understanding the structure of 
universally stable equilibria for arbitrary influence graphs, as we discuss later. For brevity, if x is 
a mass vector such that the distance between any two active nodes is at least 3, we will say that x 
is 3-separated. 

The first proposition proves the "if" direction of Theorem l4.2[ Its proof follows from our analysis 
in Section [2l 

Proposition 4.3 Under the local model with a > 1, any 3-separated mass vector x* is a Lyapunov- 
stable equilibrium. 

Proof: X* is an equilibrium by Proposition 12.11 since its active nodes form an independent set. 
By Proposition 12.31 an equilibrium whose active nodes form an independent set is stable if the 
interaction mass of each active node is greater than the interaction mass of each of its neighbors. 
For a > 1, this property holds when each inactive node has at most one active neighbor. In turn, 
this holds if and only if the distance between every two active nodes is at least 3; thus, all such 
equilibria are stable for every a > 1. D 

The "only if" direction of Theorem 14. 21 is more complicated to prove. First, we show that if the 
active nodes of an equilibrium do form an independent set, then being 3-separated is necessary to 
ensure universal stability. 

Proposition 4.4 Let x* be a universally stable equilibrium, and assume that the active nodes under 
X* form an independent set. Then, x* is 3-separated. 

Proof: For the sake of contradiction, suppose that x* is not 3-separated. Then there exists an 
inactive node u with at least two active neighbors. Let Au be the set of all active neighbors of node 
u. We use the following notation: 

s = YlveA^K'^ 77 = min^gA„<; ^ = max^gA„ 2;^ • 
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Define a = 1 + ry^. We will show that x* is not Lyapunov-stable for this a. 

To prove instability, consider a perturbation x which coincides with x* on all nodes not in 
{u\ U Aui and satisfies 

Xv < x*^ for all V ^ Au 

Xu = S for some 6 £ (0, s — /i — r/^) (7) 

(Note that s — fi — rj^ > because A^ consists of at least two nodes.) 
Under such an x, the interaction mass of node u is 

Mu = axu + y^^ Xy = r]"^ ■ Xu + s, 

v&Au 

while the interaction mass of any node v € A^ is 

My = OiXy + x„ < ax* + x„ < (1 + r/^)^ + Xu- 

Since Xu < s — ^ — rf , we have that My < (1 + rf)^! + s — ^ — rf < s, and hence, AAu > M.y. 
Thus, under this perturbation, mass starts flowing from all nodes v € A^ to n, and this continues 
until Xu > s — n — rf. Consequently, the system cannot reach any equilibrium with Xu < s — ix — rf] 
in particular, it cannot reach any equilibrium with |jx* — x|| < s — fi — rf . Since this holds for 
arbitrarily small 5, and s — ^ — rf >Oisa constant independent of 5, we conclude that x* is not 
Lyapunov-stable for this a. D 

With Proposition 23] in place, all that remains to complete the proof of the "only if" direction of 
Theorem l4.2l is to ensure that the active nodes in any universally stable equilibrium x* of a bipartite 
graph form an independent set, i.e., that each connected component C of Ta,ct{x*) consists of a 
single node. This is implied by Lemma [4.5l which shows in general that if x* is a universally stable 
equilibrium, then all the non-trivial connected components of the subgraph of its active nodes are 
not bipartite graphs. This completes the proof of Theorem 14.21 ^-s ^"^Y connected subgraph of a 
bipartite graph is a bipartite graph itself. 

An additional benefit of Lemma 14.51 is that it applies to arbitrary influence graphs, and signif- 
icantly limits the topologies a connected component of Ts,ct,{x*) can have for a universally stable 
equilibrium x* . To state this lemma in the most general form, we define a class of regular graphs 
which in particular subsumes all bipartite graphs, all cliques, and all cycles whose length is a mul- 
tiple of 3. We say that a d-regular graph is locally balanced if its vertices can be partitioned into k 
disjoint sets Fi, V2, . . . , Vfc such that each vertex v £ Vi has exactly d/{k — 1) edges to each of the 
sets Vj,j 7^ i. 

Lemma 4.5 Let x* be a universally stable equilibrium and C a non-trivial connected component 
of its active subgraph Ta.ct{x*). Then: 

(a) C is a regular graph, and x* is uniform on C (i.e., x* = x* for all u,v £ C). 

(b) C is not a bipartite graph, and, more generally, C is not locally balanced. 

Proof: We begin by proving part (a). Let li, t; € C be a pair of adjacent nodes. The equilibrium 
conditions for a = 2 imply that 2x* -|- Ylws^T ^w — 2a^u + X^iogt ^wi ^^^^ the ones for a = 3 that 
'^^v + StuGT ^*w — '^^u + Ylw&T ^*w Subtracting the first equation from the second shows that 
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xl = x* . Because C is a connected component, applying this argument along all edges in C proves 
that all nodes in C must have the same mass fi. 

The interaction mass of node v with a = 2 is therefore A4y = /i • (|T^| + 1). Considering again 
a pair u,v of adjacent nodes, the equilibrium condition A4u = My implies that \Tu\ = \Ty\. Again 
by connectivity of C, this implies that all nodes in C have the same degree, so C is regular. 

Next, we prove part (b). Because x* is universally stable, part (a) implies that C is d-regular 
for some d > 1, and x* = /i (for some //) for all u ^ C . Assume for contradiction that C is locally 
balanced, and let Vi, . . . , V^ be the k partitions of C. Because T\Vi U Vj] is a d/{k — l)-regular 
bipartite graph for each pair i ^ j, all partitions Vi must have the same size s = \C\/k. 

Set a = d + 1, and let 5 > be arbitrary. Consider perturbed vectors of the following form: 
Xt, = X* + i • (5 for every v ^ Vi and x„ = x„ — grf}_i\ ■ S for every u ^ Vi. (That is, a total mass of 
6 is removed uniformly from nodes not in Vi, and added uniformly over the nodes in Vi.) 

In moving from x* to x, the interaction mass of each node i 
0, while the interaction mass of each node u ^Vi changes by 



In moving from x* to x, the interaction mass of each node w € Vi changes by a---6—d- ^(u_i) '^ > 



1 d{k-2) 1 d 1 

"' s{k-l)' k-1 ' s{k-l)' ^AT^'s' 

d{k - 2) \ 1 



(_„,,_^^,.).^^..<„ 



s{k-l) 

Thus, for any such vector x{t) = x, all flows are directed from nodes not in Vi to nodes in Vi. 
Furthermore, by symmetry of the original vector x* and the perturbation, the new mass vector 
x{t + 1) will be of the same form, for a different 6' > 6. Thus, the same argument will inductively 
apply in every time step. Hence, the direction of flows never changes, and Lemma 12.21 guarantees 
that the system converges. Since the interaction mass of all nodes in Vi is only increasing, and 
the interaction mass of all nodes not in Vi is only decreasing, the only equilibrium y the system 
can converge to is one in which all nodes outside of Vi have zero mass. In particular, this means 
that even starting from ||x* — x|| = 5' (which would correspond to using 6 = d'/2 in our analysis), 
||x* — y\\ is not bounded in terms of 5', so x* is not stable. D 

Lemma 14.51 considerably narrows down the set of equilibria for which the question of whether 
or not they are universally stable remains open. 

More specifically, it only remains to consider mass vectors x in which there is a non-trivial 
connected C of Ts^ctix) such that C is a d-regular graph (for some d > 1), is not a locally balanced 
graph (in particular not a bipartite graph), and for every u,v G C, Xy = Xu- We conjecture that 
such mass vectors are not universally stable; it would then follow that in any universally stable 
equilibrium, all components have size 1, and hence by Proposition 14.41 the active nodes would be 
at mutual distance 3. Accordingly, we formulate the following: 

Conjecture 4.6 A mass vector is a universally stable equilibrium if and only if its active nodes 
are at pairwise distance at least 3 in the influence graph. 

4.2 Characterization of Lyapunov Stable Equilibria for a = 1 

For a > 1, we have shown that any equilibrium whose active nodes form an independent set of 
pairwise node distance at least 3 is stable. Perhaps surprisingly, this ceases to be true for a = 1. 
Indeed, on the path of length 4, the equilibrium x* = (^, 0, 0, ^) is not stable. 

We can see this instability as follows. Consider vectors of the form x^^' = (-^ — S,6,S, ^ — 5). 
Under x^ ' , for any 5 € (0, ^), the interaction mass of nodes 1 and 4 is strictly smaller than the 
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interaction mass of nodes 2 and 3 (whose interaction masses are equal) . This impHes that no vector 
x' ' can be an equihbrium for 6 G (0, 2), and that flow will always be directed from nodes 1 and 4 to 
nodes 2 and 3. Furthermore, the flow from node 1 to node 2 is equal to the flow from node 4 to node 
3, implying that at the next time step, the mass vector will be of the form x^ ' with 6' > 6. As we 
know by Theorem IB . 1 1 that the 4-path always converges, the system converges to some mass vector 
y = x^ ' such that 6* > 0. □ Since the update rule is continuous, this y must be an equilibrium. 
We have proved that the only such equilibrium is the one with 6* = 2- Thus, starting from the 
perturbation x^ ' of x*, the system can only converge to a state y in which yi =2/4 = 0. 

While a pairwise distance of 3 between active nodes is not enough to guarantee stability, a 
pairwise distance of 4 is sufficient. 

Theorem 4.7 Let x* he a mass vector whose active nodes have pairwise distance at least 4- Then, 
X* is a stable equilibrium for a = 1. 

Proof: The proof is much more involved than the proof of Proposition 12.31 for the following 
reason: even for arbitrarily small perturbations to x*, it is possible that inactive neighbors v of an 
active node u have higher interaction mass; thus, the conditions of Proposition 12.31 do not apply, 
and in fact, u could lose mass over time. However, we will be able to show that the total mass u 
loses, starting from a perturbation of magnitude at most S, is bounded by a function g{6) — )• as 

Let A be the set of all x*-active nodes, and let x(0) be a perturbation of x* with \\x* — x(0)||i < 
5 < i • miuug^x*. We will show below that for each node u ^ A, and all times t, we have that 
\xu{t) — x*! < 25. Because 

^\x^{t)-xl\ = ^Xy{t) = ^{xl-xu{t)) < ^\xu{t)-xl\, 

v^A v^A u&A u&A 

we obtain that 

||x^(t) -x*||i = ^|xu(i)-x*| < 2^|xu(t)-x*| < And -^ Q as J ^ 0. 

u&V ueA 

It remains to prove the inequality [x„(t) — x*| < 26 for all nodes u a A and times t. Define 
W = V \ {AU UugA Tu) to be the set of all nodes at distance at least 2 from all active nodes. 
We will prove the inequality by showing that any flow from u to its neighbors v can be "charged" 
against flow from W to v. More formally, we will simultaneously prove the following invariants for 
all times t and any set U '^ A: 



wgw wew 

Y ^^(*) > XI ^«(^) ~ XI (^«'(0) - Xu,(t)) for any set U C A. (8b) 

ueu ueu wew 

Let t be an arbitrary time step, and assume that the invariants hold at time t. We will establish 
that they still hold at time t + 1. First, we notice some useful consequences of the invariants, 
including the desired fact that |x„(i) — x'^l < 25. 

From Inequality dHa]), we get that XliueVK ^«'(*) ^ J2wew ^"-(O) ^ J2v(^A ^^{0) < 5. Substituting 
this bound into Inequality (|8b[) with U = {u}, and using that x«(0) > x* — (5 gives us that 



*This also follows directly from our argument with x^ '. 
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Xuit) > x"^ — 25. Similarly, using Inequality ()8bp with U = A\ {u} gives us an upper bound of 
Xuit) ^ 2;* + 26. So we have shown that \xu{t) — x* | < 25. It remains to establish the invariants 
for time t + 1. 

Let {w,v),w € W,v ^ M^ be an arbitrary edge, and u the unique active neighbor of v in 
T. The flow on the edge {w,v) is U^v{t) = p- ''"'^*^ '"^^^i^ff^^^jft)^" ^*" • We have just seen that 
X* — 2(5 < Xu{t) < X* + 2(5, so we can also bound Mvit) > x* — 26. Applying Inequality (|8bp with 
[/ = A\{u} also gives us an upper bound of 7Wi,(t) < 1 — X^u'gyi m'3^m^m'(0 < x*+2(5. Furthermore, 
using the definition of W and Inequality (\8h\i for U = A, 



MUt) < Y^x.it) < ^x,(0)+ ^ (x^(0) - x^(t)) < 26. 

Substituting these bounds, we get that 

f f,. ^ Xyj{t)x^{t){xl-4:6) 

f-^^^'^^P 26{a:l + 2^) " ^'^ 

By definition of 6, this quantity is always non-negative. In particular, this means that flow goes 
from w to v; since the edge {w,v) was arbitrary, we have established the invariant ()8ap . 
Next, fix an arbitrary node pair u (^ A,v € Tu. The flow from « to ?; is 

r (.. _ Xu{t) x^{t){Mv{t) - Mujt)) Xujt) x^{t)J2^^wnn ^Ut) 

^^^^^')-P- Mu{t)M.{t) -P- (x„(t))2 

= P 7^ Yl Xy,{t)x^{t) < P- ^ _ Y^ X^{t)x^{t). 

On the other hand, summing the bound ([9]) over all nodes w GW (IT^, we get that 

Because 6 < x*/8, we get that ^,^25 — 2(5fa*~+25) ' ^° ^^^ ^°^ from n to u is at most the total 
flow from all w £ W DTv to v. For any set U C A, summing this inequality over all n E [/ (and 
noticing that we never double-count the same edge) now establishes Invariant ()8bp at time t-|- 1. D 

5 Conclusions and Open Questions 

In this paper, we presented a novel model of cultural dynamics that captures the essential aspect of 
several previously studied models: the interplay between selection and influence. We concentrated 
on two instances of this model. In the basic version, the global model, each person selects another 
person from the entire population to interact with. In the local model, a person selects an interaction 
partner from a subset of the population consisting of similar people. 

We provided a nearly complete treatment of the global model, showing that the system always 
converges from any initial mass vector, and providing a complete characterization of stable equi- 
libria. An open question is to predict the equilibrium to which the system converges starting from 
a given initial mass vector. We suspect that with probability 1 over possible starting states, the 
system converges to an equilibrium in which the active nodes form an independent set. 

The local model involves, at its heart, a dynamical system on the population fractions that is 
complicated even for small numbers of variables. As such, it raises many interesting and challenging 
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questions, and we have made progress on some of these. In particular, we know that on paths of 
length 3 (for a > 1) and at most 5 (for a = 1), the system converges from any starting state. 
However, it is open whether convergence occurs for all graphs. On the stability frontier, for a > 1, 
we conjecture that the only stable equilibria are those in which the active nodes have pairwise 
distance at least 3. We showed that such equilibria are indeed stable, and that a number of other 
equilibria are unstable — including ones in which the active nodes form any other independent set, 
or ones in which they form a locally balanced graph (a class that includes bipartite graphs). Finally, 
we would like to raise an even more challenging question: does the dynamical system defined by 
the general model always converge? 
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A Convergence on a 3-node path (Proof of Theorem 4.1) 

Observe that the flow between two types u and v does not necessarily go in the same direction over 
all time steps, but instead may change its direction. To keep track of the changes in direction, we 
define a configuration of the system to be a labeling of all edges (u, v) in T by the direction along 
which flow is traveling (i.e., whether it travels from u to v ov from v to u). A configuration which 
cannot transition to any other configuration is called a sink configuration. Sink configurations are 
important because they guarantee convergence by Lemma |2.2[ 

In the case of a 3-node path and a > 2, there are four possible configurations. We study transi- 
tions among the configurations as the system evolves over time; we show that each configuration is 
either a sink configuration, or it has the property that any change in the direction of an edge leads 
to a sink configuration. This ensures that there can be at most one change in the direction of fiow 
as the system evolves. 

For a 3-node path and a < 2, convergence will follow as a special case of the more general 
Theorem \A.2\ which establishes convergence for all star graphs when a < 2. For a 3-node path 
and a > 2, we prove the following lemma (which, jointly with Lemma [221 implies Theorem 14. II for 
a > 2). 

Lemma A.l Consider the local model with a > 2 such that the influence graph T is a 3-node path. 
Then, there is a time to such that for any t > to, no flow changes its direction. 

Proof: Let the nodes of the path be (1,2,3), in order. Consider an arbitrary time t. Consider 
the change in the interaction mass of node 1, for example, from step t to t + 1. Recall that 
A4i{t) = axi{t) + X2(t). Node I's interaction mass is decreased by a/i->2(i) from flow leaving node 
1 to node 2, and increased by /i_j.2(t) + fz^2{'t) from fiow entering node 2. By applying the same 
reasoning to nodes 2 and 3, we get: 

Mi{t + 1) = Mi{t) + h^2{t) - (a - l)/i^2(t), 

M2{t + 1) = M2{t) + (a - l)(/l-,2(t) + /3^2(t)), (10) 

M^{t + 1) = M^{t) + /i^2(t) - (a - l)/3^2(t). 

Let Xi = Xi(t),A^i = A4i{t), fi^j = fi^j{t) for i,j = 1,2,3. We will distinguish three cases based 
on the relative sizes of TWi, A^2 5 -^3- 

1. If 7^2 ^ Ml and A^2 ^ J^3, then both fi-^2 and /3->2 are non-negative. According to 
Equation (jlOp . A^2 increases by at least as much as both Mi and A^3, so the same inequality 
will hold in the next step (and thus inductively forever). Thus, we have reached a sink 
configuration. 



2. If M2 < Ml and M2 < M3, then both /i^2 and /3_s.2 are negative. By Equation (jTO]) . 
A^2 decreases by at least as much as both A^i and Ais, so again, the inequalities will hold 
forever, and we have reached a sink configuration. 

3. The remaining case is that M2 < Mi and M2 > M3. (The case M2 < A^3,A^2 > Mi is 
symmetric.) Here, M3 decreases. Mi increases, and M2 may increase or decrease. If the 
relative order of Mi,M2,M3 stays the same for all steps after t, then we have reached a 
sink configuration. Otherwise, at some time t' > t, we must reach either a configuration with 
M2{t') < Miit'),M2{t') < M^it') or with M2{t') > Mi{t'),M2{t') > Msit'). Either of 
those configurations is a sink configuration by the preceding two cases. 

18 



In summary, each configuration is either a sink configuration, or will reach a sink configuration at 
the next transition to a different order of interaction masses. D 

Next we prove that for a < 2 the process converges on every star graph (and in particular on 
the 3-path). 

Theorem A. 2 Under the local model with a < 2, if the influence graph is a star graph, then the 
system converges from any starting state. 

In the remainder of this section we prove Theorem IA.2I More specifically, we show that even- 
tually the system enters a sink configuration. 

Lemma A. 3 Consider the local model such that the influence graph T is a star graph. Then, at 
any time, the number of edges with flow directed away from the center is at most [aj . 

Proof: Denote the central node by u. Fix some time t, and let R be the set of all peripheral nodes 
V such that the flow on the edge (n, v) is directed from u to v. Because the flow is directed towards 
V, My > A4u for all v G R. Rearranging this inequality gives us that {a — l){xv — Xu) > J2wjtu v ^«" 

which implies in particular that Xy > — "^"'^ *" ■ Summing over all v €z R now implies that 

^^-^ ^-^ a — 1 a — 1 

veR v£R 

Thus we have that \R\ < [aJ. D 

Lemma A. 4 Consider the local model such that the influence graph T is a star graph, and with 
a < 2. Then, a configuration in which a flow on exactly one edge is directed away from the center 
node is a sink configuration. 

Proof: Let u be the center node. Suppose that at time t, the system is in a configuration in 
which the flow on exactly one edge {u,v) is directed away from the center; so M.u{t) < ^Ay{t). 

We will prove that Aiuit + 1) < Mv{t + 1). By Lemma [A. 31 there can be at most one edge on 
which the flow is directed away from u, and (u, v) is such an edge; hence, the configuration at time 
t + 1 is the same as at time t. To see that Aiui't + 1) < A4y{t + 1), observe that the new interaction 
masses are 

Muit + 1) = Muit) + (a - 1) Y, fu>^u{t) - (a - l)fu^v{t) 

My{t + l)=My{t) + {a-l)fu^v{t)+ Y^ fu,^u{t). 

Now, since ^Auit) < A4y{t), their difference is 

Mu{t + 1) - X,(t + 1) < (a - 2) Y fv^^uit) - 2{a - l)/„^,(t). 

Because a < 2, the right-hand side is negative, and we have proved that M.u{t + 1) < M.v{t + 1). 

D 
Theorem lA. 21 now follows from Lemmas l2.2llA.3] and [A4l as follows. If the system ever enters a 
configuration in which exactly one edge has flow directed away from the center, then by Lemma [A.41 
it subsequently stays in this configuration forever, so by Lemma 12.21 the system converges. By 
Lemma IA.31 the only other alternative is that the system is always in the configuration with all 
edges directed inwards; then, again, it converges by Lemma [221 
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B Local model with a = 1: Convergence on a Path 

Assume that the active subgraph is an n-node path with nodes (1,2,... ,n). The endpoints of 
the path, nodes 1 and n, always have interaction masses no larger than their neighbors (nodes 
2,n — 1), implying that their masses xi(t),x„(t) are monotonically non-increasing. This implies 
convergence of xi{t), x„(t) as t — )■ oo. In the following proposition, we will exploit the convergence 
at the endpoints to show that X2{t) and Xn-i{t) must also converge. For a path of length at most 
5, this implies convergence of the vector x to an equilibrium, as the total mass stays constant. 
Our technique does not apply beyond length 5; we do not know of a direct way to generalize the 
argument inductively to paths of arbitrary lengths. 

Theorem B.l Consider the local model with a = 1. If the influence graph is a path of n < 5 
nodes, then the system converges. 

Proof: We already argued above that xi(t) and a;„(t) converge. If the path has 3 nodes, 
then X2{t) = 1 — xi{t) — x^^t) (by mass conservation), so X2{t) converges as well. So assume 
that n € {4,5}. Below, we show that X2{t) converges as well; a symmetric argument applies 
to Xn-i{t). If the path has 4 nodes, we are done at this point. If the path has 5 nodes, then 
X3{t) = 1 — xi{t) — X2{t) — X4{t) — X5(t) must Converge as well. Thus, x{t) converges in all cases. 
To prove that X2{t) converges, we distinguish two cases, based on yi = limt^ooXi{t). 

1. If yi = 0, there are two subcases. If xi{t) > X2{t) for all t, then clearly, X2{t) — )■ as well. 
Otherwise, there exists a to with 2:2(^0) > xi{tQ). By the definition of the model, specialized 
to the local model and a = 1, we obtain that for any t. 

Then, clearly, xi(t) < X2{t) implies xi{t + 1) < X2{t + 1). In particular, this means that 
X2{t) > xi{t) for all t > to- In turn, this inequality is used in the last step of the following 
derivation: 

m»(A..W,0)<p.^;||^..,(t) 

Xl{t) X2{t) X^{t) Ml{t) 

~^' Mi{t)M2{t) ' Msit) 

Xi{t)+X2{t) 



/l^2(t) 



X2{t) + X3{t) + X4{t) 
<2/l^2(t). 

Thus, the total amount of flow entering node 2 after time to is at most 3^^^ /i^2(i) ^ 
3x1 (to). The reason for the last inequality is that flow never enters node 1, so the total 
amount of flow that can leave node 1 for node 2 after to is at most the amount that was at 
node 1 at time to. 

Let F^(t) (resp., F^(t)) be the total amount of flow that has entered (resp., left) node 2 
up to time t. We have just proved that F'^{t) — F+(to) < 3xi(to). Therefore, F^{t), being 
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monotone and bounded, must converge. Because flow can only leave node 2 when it was 
already there, we get that F~{t) < X2(0) + F^{t) is also bounded, and must also converge. 
Hence, X2(i) = X2{0) + F~^ (t) — F~{t), being the difference between two convergent sequences, 
must also converge. 

2. If yi > 0, we will pursue a similar argument, but this time bounding the cumulative flow out of 
node 2 instead of into it. Because X2{t) < 1 for all times t, this means that the cumulative flow 
into node 2 must also be bounded. Then, an identical argument to the previous paragraph 
shows that 2^2 (t) = 2:2(0) + F~^(t) — F~{t) must converge. 

No flow can ever leave node 2 for node 1, so we just need to bound the flow from node 2 to 
node 3. Flow leaves node 2 for node 3 at time t if and only if Xi{t) > xi{t). Since xi(t) > yi, 
it follows that M2{t),-M3{t) > yi as well. Therefore, we can bound the non- negative flow 
from node 2 to node 3 as follows: 

If f.\ n\ ^ X2{t)x3{t) X2it)x3it) 
max(/2^3(t),0)<p.-^^-^^^^^^.X4(t) < p -^ 

Xl{t) X2{t) X3{t) 1 Xl{t) X2it) X3{t) 1 

-' 7i 7,-'- M,{t)M2{t) = ^•^^-(*)- 

Thus, the total positive flow from node 2 to node 3 is bounded above by a constant times the 
total flow from node 1 to node 2, which in turn is at most xi(0). 
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